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Abstract 

We study the so-called nonmagnetic phases (dimer and flux states) 
in the t-J model below half filling. We present a new phase diagram, 
at zero and finite temperature, that includes broad areas of phase co- 
existence (dimer-flux or flux- uniform) , in accordance with experimental 
and numerical data on the possibility of separation into hole-rich and 
hole-poor regions. We also briefly comment on some techniques used in 
the literature to discuss phase separation in the t-J model. 

PACS: 71.27, 74.20. 



Since Anderson's suggestion that strongly correlated electron models 
might be relevant to describe the physics of high-T^ superconductors, much 
work, both numerical and analytical, has been done in order to understand the 
properties of the Hubbard model in 2-dimensions, in particular by examining 
simplified models that can be obtained from it in the strong coupling limit, 
such as the Heisenberg model or the t-J model for the one-band version or 
the spin-fermion model for the three-band case. 

Still, the question of the ground state of these models has not been settled, 
at least away from half-filling. As candidate for the ground state, a variety of 
phases has been proposed, which include the long-range antiferromagnetically 
ordered (AF) state 0], the short-range antiferromagnetically ordered (RVB) 
state 1^, as well as the so called non-magnetic states, such as the dimer and 
the flux phase 0. 

In this paper, we will concentrate on the t-J model, which has been exten- 
sively studied by means of different techniques and within different approaches, 
that include standard mean field theory [6-8], slave-boson techniques [9-12] and 



numerical simulations [13, ^ 



At half-filling, i.e. at zero doping {6 = 0), when the model reduces to the 
well known spin 1/2 Heisenberg Hamiltonian, the instability of the other phases 
against the AF phase has been established beyond any reasonable doubt by 
analytical and numerical methods |T^. Within mean field theory, the lowest 
energy phase is not the AF but the dimer phase, either the columnar or the 



staggered one [|T6], which are degenerate at this level of approximation. 



Even the inclusion of both classical and quantum quadratic fluctuations does 
not increase the energy of the dimer phase above the energy of the AF state 
|T8[] . It is known that this apparent contradiction is strictly related to the 
difficulty of taking into account the constraint of one electron per site in this 
approach. Indeed, in |jl8| it has been shown that at the mean field level. 



eventually corrected by quadratic fluctuations, the constraint is implemented 
exactly for the AF phase but not for the dimer phase. Thus, in the latter case 
we are effectively dealing with a larger space of states and hence it comes of 
no surprise the fact the energy gets substantially lowered. 

For finite doping fraction (5 > 0), the major difficulty consists of developing 
good approximating techniques while dealing with the Gutzwiller (or below- 
half-filling) constraint. The experimental evidence and all the above mentioned 
theoretical studies confirm the intuitive idea that the presence of holes tends 
to disfavour a long-range antiferromagnetically ordered phase. However, there 
is still no complete agreement on the possible ground state at 5 7^ 0. Among 
the candidates for the ground state of the t-J model, the RVB and the flux 
states have been mainly considered. The former, first proposed by Anderson 
and coworkers [0 , is characterized by a superconducting ordered parameter of 
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the kind A,j = {ci^Cji — CiiCji^) , where z, j are nearest neighborhood sites, which 
gives short-range antiferromagnetic correlations. The latter corresponds to 
a nonmagnetic complex order parameter Uij = (Sa=T,j, ^la^i")' whose phase 
originates a nonzero magnetic flux threading the elementary (square) plaque- 
ttes of the lattice in a staggered way. Whenever the flux is different from or 
TT (mod 27r), this phase breaks parity and time- reversal symmetries. This has 
been the starting point of a large number of works, in which it has been argued 
that the continuum limit effective action describing the low-energy excitations 
around the flux phase could contain a Chern-Simons term [1^, which would 
impart fractional statistics to the quasiparticles, leading to the possibility of 
anyonic superconductivity ||2^ . An analysis of this claim requires more sophis- 



ticated techniques than the one used in this paper and goes beyond our present 



scope. We remark here only that, while it exists for the AF phase pl| , ^, a 
detailed study of the continuum limit around the flux phase, able to prove or 
disprove the presence of a topological term, is still missing in the literature. 
We plan to further discuss this point in some future work. 

At half flUing, the large-U limit of the Hubbard model admits an SU(2) 
gauge symmetry , which shows the equivalence between the apparently dif- 
ferent RVB and flux states. This symmetry is broken at 5 7^ 0, so that these 
two states evolve into different phases, which nevertheless stay very close in 
energies. In this paper we will not address the problem of the competition 



between the RVB and the flux phase We will concentrate our analysis 

on the nonmagnetic phases only. This is because we seem to flnd a phase 
diagram in the temperature-doping (T-S) plane richer than the one considered 
in previous works within the same kind of mean fleld approximation. If, on 
one side, our pure mean fleld data conflrm the hypothesis that the flux phase 
gives the ground state in the range of the doping fraction 6 which is relevant 
for superconductivity, on the other side, we discover that, especially at low 
temperatures, there are wide regions of phase coexistence, corresponding to 
either dimer-flux or flux-uniform phase separation, depending on the value of 
5. This seems to go in the same direction of some experimental data on high 
temperature superconducting materials P3, |25| and as weel as of some pertur- 



bative and numerical studies of the t-J model [26-30], which flnd a tendency 
towards a separation of regions with different hole concentrations. We remark 
here that the presence of phase coexistence regions might also affect a stability 
analysis of the RVB vs. the flux phase. 

Our starting point is the t-J Hamiltonian, which can be obtained as the 
strong coupling limit, \t\ <^U, of the one-band version of the Hubbard model 
for small values of the doping 6 : 

"^e// = J2 Uj{l-nia)cl^Cja{l-nja) (l) 
<ij>,a 
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<ij> ^ 

where c|q, is the creation operator of one electron with spin a =t, i in the site 
i and 5'j = | J2af3 cl^^aisCip (c? = (ai, (T2, Us) being Pauh matrices) are the spin 

operators on the site i. In addition, Jij = ^^^f^-, a = —a while the symbol (ij) 
denotes a sum over nearest-neighbor (n.n.) lattice sites only. The operator 
expressions containing riia = c\^Cia have the effect of enforcing the Gutzwiller 
projector. 

At exact half- filling, when = rii^ + = 1, the first term of (|l]) van- 
ishes identically, while the second line reduces to the spin ^ antiferromagnetic 
Heisenberg Hamiltonian. If we introduce holes in the system, so that the dop- 
ing fraction is 5 7^ 0, we can no longer neglect the first term of (|^) which 
describes a direct hopping of one electron from the site j to the (empty) n.n. 
site i. Because of the presence of the number operators rijo,, the hopping term 
is the sum of products of up to six electron creation/annihilation operators 
and is therefore very difficult to analyze. 

Thus, following Anderson and coworkers 0, we make the assumption that 
the main effect of the Gutzwiller projection is the renormalization of the hop- 
ping amplitude from its nominal value tij to 5tij. We recall that such assump- 
tion has been proven to be correct either within a slave boson approach at the 



mean field level |31| and within a variational approach making use of the 
Gutzwiller approximation ||3^. With the further assumption that tij = t, up 



to nonrelevant constant factors, the Hamiltonian (|T]) becomes: 

^ = ^'^ E E 4aCja - ^ EE 4aCjac]pCip . (2) 
<ij> a ^ {ij) al3 

We will also restrict our attention to 2-D square lattices. 

To study the phase diagram of (Q), we will follow the technique described 
in [|16[ to rewrite the partition function in the grand canonical ensemble Z = 
Tr as: 



Z = |[P^*„P^,„]exp|- J^'dr 



EE^taidr - ^^ma (3) 



+ /^^^ E E ^L^ia - ^ E E ^ia^j»^jf3^i(^ 
<ij> OL ^ <ijy a/B 

where we have introduced Grassmann fields ipia for the fermionic operators 
We can now decouple the quartic term in the exponential of (0) via 
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Hubbard- Stratonovich auxiliary fields lAij as follows: 

Z = j [Vi;l V^,^] I [VU:pU.,] exp I - ^ t/r ^ [ (4) 



<ij> 



nl 

X exp { — j dr 



i OL 



<ij> a <:ij> i3 

We will work in the static approximation: Uij{T) = Uij (constant in the 
imaginary time r), so that, going to Fourier transform with respect to r, (|^) 
can be rewritten as: 



(5) 



X exp I E E E - (3v]t>^,^M 

nn' ij ct 

where a;„ = (2n + l)7r [n G 2) are Matsubara frequencies and 

= («a;„ + /i/?)(5,/„,„, (6) 

ri|Ti*i ^ \ \tb ^Uij\bnn' ifin.n. j , . 

^ J"*^' I otherwise ' ^ ' 

We can now perform the integral over the fermionic variables to get: 

Z = |[PW*.P%]exp|-^5:W*.W,,jexp{-5,y^ (8) 

S,ff = -2Tr{log[-G-i+/3U]} , (9) 

where the factor 2 comes from spin summation and "Tr" stands for a trace 
over lattice sites and frequencies. 

Following we consider solutions of the saddle point equations, = 
0, that admit a symmetry for translations along the diagonal of elementary 
plaquette of the square lattice. In this case, the matrix U depends only on four 
independent link variables Uj as shown in figure |l|(a). Under this assumption, 
the corresponding Brillouin zone gets halved. The Reduced Brillouin Zone 
(RBZ) is given by the shaded area of figure |l|(b). 

In momentum space, the matrix U can be easily diagonalized. It has eigen- 
values: 

Ek = ± lAfel (10) 




Figure 1: 

(a) The periodicity of the four independent link variables Uj, j = 1, ■ ■ ■ ,4:. 
(b) The hatched area shows the Reduced Brillouin Zone. 



with Xj =t6 + Uj and k e RBZ, so that: 

keRBZ n 

(11) 

The sum over the Matsubara frequencies can now be performed, yielding the 
following expression for the partition function in the static approximation: 

Z = j [VU*VUj\e-'^^ (12) 
Q = -jY.W-^ E [log(l + e'^('^-l^^l))+log(l + e''('^+l^^l)) . 

j=l keRBZ 

The stationary points of are thus given by the saddle-point equations: 



J 


E 

keRBZ 


-ik^a ■^k 


1 


1 




, (13) 


iV 






1 + g-^^/'+l^fcl) 




J 


E e- 

keRBZ 


A* 


1 


1 




, (14) 


N 








1 + e-^(M-|Afc|)_ 


J 


E 

keRBZ 




1 


1 




, (15) 


N 




1 + e-^(M+|Afc|) 


1 + g-^(M-|Afc|)_ 
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* AT ^ 



^ k&RBZ l-^^l 



1 1 



1 _^ e-/3(M+|Afc|) 1 _^ e-/3(M-|Afc[) 



(16) 



which must be supplemented by the equation determining the chemical poten- 
tial ii: 



N{l-8)=2 Y. 

k&RBZ 



1 

+ 



1 _^ e-/3(Ai-|Afe|) ]_ _^ g-/3(/i+|Afc|) 



(17) 



We have already recalled that, at half filling, the Hamiltonian (§) admits an 
SU{2) gauge symmetry. As a special case, (0) is invariant under a local U{1) 
transformation mapping the operator Cja into e^^^Cja- Such transformation 
changes the phases of the link variables Ujk = \lAjk\e^^^'' by 

9,jk — > 9jk + 6k — 9j , (18) 

but keeps the sum of the phases around an elementary plaquette (of vertices 
j, k, I, m) invariant: = 6jk + Oki + Oim + Omj = cost. Indeed, A9 is a gauge 
invariant , hence observable, quantity and can be thought of as the fiux as- 
sociated to a magnetic field threading the plaquette. In addition, equations 



(p!^-|T7|) with 5 = are left invariant and hence the grand canonical potential 
Q will admit a set of degenerate minima, parametrized by U{1) and all corre- 
sponding to the same A^. For S ^ 0, this symmetry is explicitly broken by the 
hopping term and, in general, only one specific choice of the phases 6jk, which 
are no longer pure gauge degrees of freedom, will correspond to a minimum of 

n. 

We have chosen to study equations (|T^-|T^) numerically for the case of 
t/J= 1, a value that well approximates the experimental values for the hop- 
ping amplitude (t = leV) and the Hubbard repulsion {U = 5eV). We have 
found the following solutions yielding (local) minima of the grand canonical 
potential Q: 



1. Uniform phase. 

This is characterized by the choice: 

Ui=U2=U^=Ua=U eWi^ . (19) 

This solution exists at all temperature and is the only minimum of Vt at 
high temperature and/or high doping. At all value of the doping 5, the 
total fiux per plaquette is /S.6 = 0. 

2. Dimer phase. 

It corresponds to solutions of the form: 

Ui^U2=U^=Ua, Uj G (D (20) 
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and with 

\Ui\ > IW2I • (21) 

At (5 = 0, the ?7(l)-gauge symmetry mentioned above allows to choose the 
Uj real. In addition, in this case, U2 = 0. For finite 6, gauge invariance 
is broken and the parameters Uj are in general complex. This is indeed 
the case for a range of temperature < T < Ti = O.IJ, in which the 
parameters stay real only up to some value (5ci(T) and then acquire a 
phase increasing with 6. For T > Ti, on the contrary, the parameters 
are always real. In addition, for 6 > Sc^iT), (^Dp ceases to be a solution 
of the saddle point equations. This behaviour is shown in figure 




Figure 2: 

The phase of the order parameter Ui in the dimer phase 
as a function of the doping S at different temperatures. 
The critical dopings Sci and Sd are shown for T = 0.02. 



3. Kite phase. 

This solution exists only in a very small region of the parameter space 
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and corresponds to the choice: 

Ui=U2^U^ = Ui ,Uj eWi . (22) 

The parameters and W4 are zero at 5 = and then increase while 5 
increases. 

4. Flux phase. 

This phase corresponds to a sort of complex uniform phase and is given 
by the choice: 

Ui=U2=U^=Ui= Ue''*' , W e 1R+ and ^ . (23) 

Because of the U{1) invariance, at 5 = 0, the phases of the four param- 
eters do not need to be the same, as long as the total magnetic flux per 
plaquette is equal to 40. In figure 3, is shown as a function of 5 at 
different values of the temperature. For any T < T^^\ is | at exact 
half-filling and then decreases with continuity going to zero at some crit- 
ical value of delta 5c{T). We will discuss the behaviour of the curve for 
T = 0.02 shortly below. 



After having calculated and compared the free energy F = VL + fiN{l — 6) 
corresponding to the different mean field solutions described above, we have 
obtained the phase diagram shown in figure 4. We remark that, contrary 
to what happens in going from the dimer to the kite or the flux phase, the 
transition from the flux to the uniform phase is continuous, the phases of the 
order parameters Uj going smoothly to zero as 6 increases. 



We observe that, for T = 0, we recover the results of p. Indeed, since we 
have renormalized the hopping coefficient t by a factor of 6, our data have to be 
compared with the ones that can be read moving along the diagonal line t/J = 
6 of the zero temperature phase diagram given in [^. Also, we remark that the 
flux phase has been extensively studied [^, |Tl| by analyzing the Hamiltonian 
(0) within a slave boson approach. Apparently, some different results are 
obtained: the authors of references |12|, |11| find a fiux phase corresponding 
to a magnetic fiux per plaquette which for small 6 is blocked to the value 
vr, decreasing towards zero only after a finite value of the doping. It is not 
difficult to show, however, that the Hubbard- Stratonovich fields U-j used in 
those papers differs from the Uij we have used, the relationship being: 
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Figure 3: 

The phase cj) of the order parameter in the flux phase 
as a function of the doping S at different temperatures. 



where bi is the creation (annihilation) operator of the slave boson and the 
last equality is valid in the mean field approximation. We have checked that 
our numerical data for the Uij in the flux phase, once corrected according to 
formula (|2^), yield a phase for the U-j that does indeed reproduce the behaviour 



described in |T^, |TT[. Thus, as it should be, the two approaches give the same 
mean field results. 

Our mean field calculations indicate the existence of a rather extended 
region of stability for the flux phase below half-filling, which coincides roughly 
with the region relevant for superconductivity. We notice also that, between 
dimer and flux as well as between flux and uniform phases, there is a rather 
strange receding of the phase boundary line. This is the cause of the behaviour 
of the T = 0.02 curve in figure 3. 

The phase diagram shown in figure 4 has been constructed without taking 
into account the possibility of phase separation. However, this possibility 
cannot be overlooked. It is indeed a long-standing and still unresolved problem 
concerning the Hubbard model and its parent Hamiltonians, such as the t-J. 
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Figure 4: 

The phase diagram of the t-J model at pure mean filed level. 



Let US recall that experiments on different copper oxyde compounds have found 
a spatial separation between AFM hole-poor and superconducting hole-rich 
regions, happening either at a microscopic level p4| or at a macroscopic level, 
such as in ordered striped phases p5|. From a theoretical point of view, much 



work has been done to demonstrate that this phenomenon can occur for the 
Hubbard model in the strong coupling limit. In particular phase separation for 
the t-J model has been discussed by means of a high-temperature expansion 

by using variational arguments 



by Puttika et al. p6[ and by Emery et al. 
and exact diagonalization numerical studies. Its connections with the onset of 
high-temperature superconductivity has been considered in |Q. Also, phase 
coexistence has been examined for other phases that might be relevant in the 

H 



t-J model, such as the spiral ones 



or the ferromagnetic one [^, 0, in 
the latter case for its relevance to the Nagaoka problem. Anticipating our 
results, we find, for nonmagnetic phases as well, an instability towards phase 
coexistence and separation between hole-rich and hole-poor regions. Since 
we work in a mean field approximation, we are however not able to say how 
spatially distributed the two different regions can appear. 
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We have therefore inquired about the stabihty of our mean field solutions 
towards phase separation by studying: 

1) The fluctuation matrix for the grand canonical potential fl around the 
saddle point solutions, since the mean field solution is unstable whenever the 
fluctuation matrix has one or more negative eigenvalues. 

2) The behaviour of the chemical potential function of doping 6, since 
an instability towards phase separation is sig naled by |^ > |2§. In fact, if 

> then the isothermal compressibility — ^'^^ [n = 1 — 5) is negative 
and the system is clearly unstable. The stable state is given by the coexistence 
between two phases having different values of doping and the same chemical 
potential. These values of doping are given by the application of the Maxwell 
construction to the diagram . 



3) The convexity of the free energy density / = F/V as a function of 5, since an 
instability towards phase separation is indicated also by a region of concavity of 
the free energy /. It is easy to see that, when working within an approximation 
scheme that respects the thermodynamical equalities, the region of concavity 
of / should coincide with that where |j > 0, since = = n^|^. In this 

case it is also possible to check that the bitangent to the curve f{6) gives an 
interpolation between two different phases which is equivalent to the Maxwell 
construction for fi{S). 

As for 1), we have considered only fluctuations around the mean field so- 
lutions which preserve the symmetry under translations along the diagonal of 
the square elementary plaquettes. We have found that the fluctuation matrix 
develops a negative eigenvalue in two regions of the phase diagram: a) on the 
left of the boundary separating the flux and the uniform phase, below a cer- 
tain temperature ~ 0.145 J; b) on the left of the boundary between dimer 
and flux phase, below a given temperature Td ~ 0.085J. We might there- 
fore expect both a region of flux-uniform and a region of dimer-flux phase 
coexistence. 

Such scenario is fully supported by the behaviour of the chemical potential 
as a function of 6. Indeed, || becomes positive exactly in the regions where 
the fluctuation matrix develops a negative eigenvalue. We will analyze in detail 
only the flux-uniform phase coexistence, the discussion for the dimer-flux case 
being completely analogous. 



The function fi{6) for the flux and uniform phases is represented in figure 
5: for 5 below the critical Sc{T) of the transition between the flux and the 
uniform, we have plotted the chemical potential of the flux phase, while for 
S > Sc{T) we have plotted the chemical potential of the uniform phase (this 
explains the cusp at 5c(T)). While for T > Tc ^ 0.145J the chemical potential 
is always decreasing with 6, when T < TcWe have |t > in a region of doping 
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Figure 5: 

The chemical potential as a function of the doping S for the flux and the uniform phases. 



which extends from a certain 6co{T) up to Sc{T). In this region the flux phase 
is, therefore, unstable and the system separates in a hole-rich part, which is in 
the uniform phase, and in a hole-poor part which is in the flux phase. 

The thermodynamical identities would require the free energy to be concave 
only in the region of the S — T diagram with Sco{T) < S < Sc{T) and T < Tc. 
But we have found that, for temperature above Tq ~ 0.12 J, the free energy 
of the flux phase is already concave at 5 = and remains concave up to the 
critical Sc{T). We believe that such thermodynamical inconsistency is entirely 
due to some shortcoming of the mean field approximation. In fact, let us note 
that the inclusion of thermal and quantum fluctuations might improve the 
consistency between the convexity of f{S) and the monotonicity of fi{S), since 
it could substantially modify the behaviour of the free energy while keeping 
the chemical potential unchanged, the latter being uniquely determined by 
the saddle point equations. In the case we are dealing with, such inclusion is, 
however, nontrivial and somewhat problematic. Indeed, at 5 = 5co{T). one of 
the eigenvalues of the fluctuation matrix becomes zero, so that the integration 
of gaussian fluctuations would give actually a divergent result. In a similar way. 
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special care has to be taken in the region of small 6: at 6 = the fluctuation 
matrix has, for the flux-phase, three zero modes corresponding to the U{1) 
gauge invariance discussed above. At exact half filling, one can get rid of these 



modes |T8[ by the standard Faddev-Popov procedure. Below half filling this 
gauge invariance is broken, but for small 6 the eigenvalues are almost zero and 
a simple gaussian integration would overestimate the fluctuations. It is clear 
that next to quadratic corrections would be required to obtain a meaningful 
result. 

We have not found any mention to inconsistencies between the behaviour 
of the free energy and that of the chemical potential in literature. In all the 
previous articles discussing phase separation in the t-J model, people have 
considered only one of the two approaches, but we believe that, at least for 
the works dealing with the mean field approximation, problems like the ones 
considered above would have to be found. We have been able to exphcitly 
check this in one example, by comparing the Maxwell contruction on the free 
energy which is given in |^ with the behaviour of the chemical potential, 
whose analytic expression is also reported. 

Let us remark also that, at temperature T > To, we find the free energy to 
be strictly concave in the range < 6 < 6c{T), which does include the zero. 
It is easy to see that this problem cannot be overcome by going to negative 
values of the doping (i.e. above half filling), since two phases with 6 of opposite 
sign cannot have the same chemical potential. Thus, in some papers p8[] , 
the "bitangent construction" on the free energy is performed by drawing the 
tangent to the curve f{6) from the point (0, f{S = 0)), forgetting the fact that 
this line is not really tangent at 5 = 0. But it is not difficult to check that this 
leads to considering the coexistence of two phases having different chemical 
potentials, which cannot therefore be in equilibrium. 



For all the reasons explained above, we have decided to base our analysis 
of phase separation on the behaviour of the chemical potential only, which, 
besides being stable against the inclusion of fluctuations, is also in agreement 
with the data on the positivity of the eigenvalues of the fluctuation matrix. As 
we have already said, we find two regions of instability, corresponding to the 
possibility of coexistence between a) the flux and the uniform phases, b) the 
dimer and the flux phases. The Maxwell construction on the chemical potential 
and the evaluation of the free energy for the mixed phases lead to a new phase 
diagram, which is substantially different from the one of figure 4. At high 
temperatures, the phase separation instabilities disappear and hence figure 4 
does not change. On the contrary, at low temperatures, the regions of dimer- 
flux and flux-uniform phase coexistence are quite broad, and for T < 0.05J 
the pure flux phase is never a minimum of the free energy, for any value of the 
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Figure 6: 

The phase diagram for the t — J model, after taking into account phase separation. 



doping. In addition, the kite phase is no more present, since it has a higher 
energy than the coexistence between the flux and the uniform phase. The 
final phase diagram, taking into account phase coexistence, is shown in figure 
6. Notice that we do not find receding phase boundaries anymore. 
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